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ABSTRACT 

We present observations of the first two supernovae discovered with the recently installed Advanced Camera for 
Surveys (ACS) on the Hubble Space Telescope {HST). The supernovae were found in Wide Field Camera images 
of the Hubble Deep Field North taken with the F775W F850LP, and G800L optical elements as pait of the ACS 
guaranteed time observation program. Spectra extracted from the ACS G800L grism exposures confirm that the 
objects are Type la supernovae (SNe la) at redshifts z=0.47 and z=0.95. Follow-up HST observations have been 
conducted with ACS in F775W and F850LP and with NICMOS in the near-infrared FllOW bandpass, yielding 
a total of 9 flux measurements in the 3 bandpasses over a period of 50 days in the observed frame. We discuss 
many of the important issues in doing accurate photometry with the ACS. We analyze the multi-band light curves 
using two different fitting methods to calibrate the supernovae luminosities and place them on the SNe la Hubble 
diagram. The resulting distances are consistent with the redshift-distance relation of the accelerating universe 
model, although evolving intergalactic grey dust remains as a less likely possibility. The relative ease with which 
these SNe la were found, confirmed, and monitored demonstrates the potential ACS holds for revolutionizing the 
field of high-redshift SNe la, and therefore of testing the accelerating universe cosmology and constraining the 
"epoch of deceleration." 
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The discovery of an apparent acceleration in the universal 
expansion (Riess et al. 1998; Perlmutter et al. 1999) using 
Type la supernovae (SNe la) has been one of the prime moti- 
vations for the recent dramatic shift in the prevailing cosmo- 
logical paradigm. The formerly favored, spatially flat, matter- 
dominated Einstein-deSitter model has given way to the new, 
observationally mandated "concordance model" possessing a 
higher expansion rate, a lower matter density consistent with 
long-standing indications from cluster and large-scale structure 
studies, and a mysterious dark energy component usually iden- 
tified with Einstein's cosmological constant A. Like the de- 
posed paradigm, the new one is also spatially flat with a total 
energy density equal to the critical density, as required by cos- 
mic microwave background (CMB) analyses (e.g., Jaffe et al. 
2001; Pryke et al. 2002), but the ratio of dark energy to matter 
at the present epoch is roughly 2: 1 . 

While a dark energy component is the favored explanation 
for the unexpected faintness of distant SNe la, other possibili- 



ties include luminosity evolution and grey dust obscuration. In 
their simplest form, these alternative models predict that SNe la 
should continue to grow fainter with respect to a fiducial uni- 
verse having constant expansion. On the other hand, dark en- 
ergy models predict that SNe la should become differentially 
brighter at redshifts z > 1 as a result of past deceleration. So 
far, only the uniquely distant SN 1997ff at z ~ 1 .7 has provided 
evidence for past deceleration (Gilliland et al. 1999; Riess et al. 
2001), and this object was never spectroscopically confirmed as 
Type la (although its host and colors indicate it was). Only fur- 
ther high-quality data on SNe la at z > 1 can confirm the cosmo- 
logical explanation and eliminate the competing astrophysical 
ones. Such data will also help in pinning down precise values of 
the cosmological parameters, which are now only weakly con- 
strained without the inclusion of external information, such as 
from the CMB. 

Fortunately, the hunt for high redshift supernovae is becom- 
ing easier. Large format CCD mosaic cameras can efficiently 
search over a much wider area of the sky and spectra of faint 
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SNe la candidates can now be obtained using 8-meter class tele- 
scopes. Another recent development expected to benefit the 
field of high-z SNe la research is the installation of the new Ad- 
vanced Camera for Surveys (ACS) (Ford et al. 1998) on the 
Hubble Space Telescope (HST). ACS has roughly twice the field 
of view, higher resolution, and nearly five times the sensitivity 
of WFPC2 at 8000 A; this combination promises to increase the 
discovery efficiency of HST by an order of magnitude. 

Here, we present observations of the first two SNe la dis- 
covered with ACS. The objects were found serendipitously in 
images of the Hubble Deep Field North (HDFN) taken by the 
Instrument DevelopmentTeam (IDT) as part of the ACS guar- 
anteed time observation (GTO) program. The following section 
describes our discovery and follow-up observations, extraction 
of the grism spectra, and photometric measurements. Section 3 
presents the SNe la Ught curves, the inferred distance moduli, 
and the updated Hubble diagram. Section 4 then discusses the 
implications of our results for cosmology and the prospects of 
ACS for the high-z SNe la search field. 

2. OBSERVATIONS AND DATA REDUCTIONS 
2.1. Imaging Data 

The Hubble Deep Field North (HDFN) was observed with 
the ACS Wide Field Camera (WFC) as part of the GTO pro- 
gram, proposal number 9301. The observations were conducted 
on 2002 May 11 in the F775W (SDSS and G800L (grism) 
bandpasses for two and three orbits, respectively, and on 2002 
May 21 in F850LF (SDSS z) for three orbits. A complete de- 
scription of the GTO program will be given in a future paper; 
here we concentrate on our measurements of the two super- 
novae that were found in these data. 

The images were initially processed through the STScI 
CALACS pipeline (Hack 1999), which performs the bias and 
dark subtraction, flat-fielding, and conversion of the counts to 
electrons. Further reduction was done through the "Apsis" GTO 
data pipeline, described in detail by Blakeslee et al. (2002). 
Apsis measures offsets, rejects cosmic rays (CRs) and detec- 
tor defects, and combines the images to a single geometri- 
cally corrected one using the drizzle method (Fruchter & Hook 
2002). We used the default "square" kernel (bilinear interpola- 
tion) with an output scale of 0."05/pix for the drizzling. The 
resulting full width at half maximum (FWHM) of the point 
spread function (PSF) was about 0."104 in F775W and Of 110 
in F850LP. While alternate kernels or output scales could re- 
duce the PSF FWHM, this approach allows us to apply directly 
the aperture corrections from the ACS photometric calibration 
program (Sirianni et al. 2003). 

Inspection of the reduced images revealed two bright point 
sources that were not present in the original 1995 HDFN im- 
ages. Figures 1 and 2 show subsections of the F775W and 
F850LP ACS images centered on these objects, as compared to 
the same regions taken from the 1995 WFPC2/F814W images. 
The first object, designated SN2002dc (Magee et al. 2002), is 
located at R.A. = 12''36'"49?84, Dec. = -f62°13'13^'0 (J2000), 
where our coordinates are on the system of the Fernandez-Soto 
et al. (1999) catalogue. SN2002dc is 0."84 from the center 
of its apparent host galaxy, which is number 2-264.1 in the 
Williams et al. (1996) HDFN Catalog and has a measured red- 
shift z = 0.475 (Cohen et al. 1996). 

The second object, SN2002dd (Tsvetanov et al. 2002), was 
originaUy identified by its spectrum; it has coordinates R.A. = 
12'>36'"55.*36, Dec. =-f62°12'46:'1 (J2000). It is located l."64 



from the center of the HDFN Catalog galaxy 3-486.0; how- 
ever, this is not the host since it has a redshift z = 0.79 (Co- 
hen et al. 1996), whereas the SN2002dd spectrum (see below) 
places it at z » 0.95 (note that the 11% higher redshift from 
the discovery telegram was based on the rough ground-based 
calibration of the ACS grism). Interestingly, the diffuse galaxy 
3-443.0 has z = 0.950 (Cohen et al. 2000) and is centered 3'.'6 
away (projected separation < 20 /;"' kpc). In addition, there is a 
fainter, uncatalogued object 0."87 (4.7 /z"' kpc) from SN2002dd 
projected between it and 3-486.0 (see Figure 2). This object's 
grism spectrum is consistent with a starburst galaxy at z=0.95 
(see the following section), and there is some indication of 
faint luminosity (possibly a spiral arm) arching from it towards 
SN2002dd. We consider this the most likely host galaxy. 

Because the ACS grism spectra indicated that these objects 
were both SNe la near maximum luminosity, we reobserved this 
field on 1 3 June under HST GO proposal #9352 for one ACS or- 
bit split between F775W and F850LP and one orbit with Cam- 
era 2 of the recently revived NICMOS instrument in the Fl lOW 
bandpass (SN2002dd only). Two additional single-orbit ACS 
observations done under GTO proposal #9301 on 24 June (split 
between F775W and F850LP), and 30 June (F850LP). These 
follow-up images were processed in the identical manner as the 
original HDFN GTO data, except that they were CR-spUt and 
CALACS performed the CR rejection. The NICMOS images 
were processed through the STScI CALNIC software and then 
combined using drizzle. 

2.2. Grism Spectra 

PreUminary discussions of the ACS G800L grism perfor- 
mance are given by Walsh et al. (2002) and Ford et al. (2002). 
The spectral dispersion for WFC grism exposures is 40 A pix"' 
on average, but varies with position by ±10%. The resolu- 
tion is typically ~90 A. Spectra of the two HDFN supemovae 
were extracted from the reduced ACSAVFC grism images us- 
ing the aXe software (Pirzkal et al. 2002), which also performs 
the wavelength and flux calibrations. The on-orbit wavelength 
calibration is accurate to < 10 A, or ~ 0.001 in redshift, and has 
been discussed in detail by PasquaU et al. (2002). We subtracted 
a polynomial fit to the sky in the cross-dispersion direction in 
the 2-d image before aXe extraction. For SN2002dc, we also 
needed to fit and remove the contribution from the host galaxy, 
so we masked the region of extraction and fit a polynomial to 
the very close background, interpolating over the masked re- 
gion. The resulting local background was very flat. 

We cross-correlated the extracted spectra against a database 
of high-signal-to-noise spectra of 27 different SNe la, most of 
which had data from multiple different dates. We determine 
the redshift from each template that gave a significant and un- 
ambiguous correlation peak, taking the rms dispersion in red- 
shift values as an estimate of the uncertainty. We also per- 
formed cross-correlations against spectral templates for other 
supemovae types (Filippenko 1997), but none of these yielded 
an acceptable fit to either of the spectra. For SN2002dc, we 
find excellent matches with normal SNe la such as SN1995D, 
SN1994ae, SN1995E, and SN1994S. These overlap more flian 
half of the rest-frame spectral range and give cross-correlation 
r values (Tonry & Davis 1979) of 16 or more. The best consis- 
tency is for an age f = 1 ± 3 days past maximum, and the best- 
fit redshift is z= 0.473 ± 0.006, with a total range from 0.45 to 
0.49. This is in excellent agreement with the z = 0.475 redshift 
for die host galaxy from Cohen et al. (1996, 2000). 
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Fig. 1. — Comparison of a small region of our ACS images centered on the position of SN2002dc with the same region from the original 1995 WFPC2/F814W 
image of the field. The displayed field size is 12"5 X 12^5 (250 ACS WFC pixels on a side) and the position angle is 140° . The separate panels show: ACS/F775W 
image (a); ACS/F850LP image (b); original WFPC2/F814W image (c); ACS F850LP image with the scaled WFPC2/F814W image of the SN2002dc host galaxy 
subtracted (d). The ACS images are displayed without any relative scahng of the counts, showing the closely comparable signal levels in the two bands. The three 
galaxies in this field with measured redshifts (Cohen et al. 2000) are labeled in the lower left (WFPC2) panel, including the z=0.475 host galaxy. SN2002dc is 
immediately apparent after subtraction of the host. 
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Fig. 2. — Comparison of our ACS and NICMOS images centered near SN2()02dd with the same region from the original 1995 WFPC2/F814W image of the field. 
The displayed field size and orientation is the same as Fig. 1. The panels show: ACS/F775 image (a); ACS/F850LP image (b); original (1995) WFPC2/F814W 
image (c); follow-up NICMOS/Fl lOW image one month after discovery. Galaxies in this field with measured redshifts (Cohen et al. 2000) are labeled in the lower 
left (WFPC2) panel. The candidate host galaxy identified in the text is labeled in the upper left panel; it has AB magnitudes F775W~25.5, F850LPrj25.3. 
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Fig. 3. — ACS grism spectra transformed to the rest frame for SN2002dc at 
z=0.47 and SN2002dd atz=0.95 are compared with ground-based spectra from 
two redshift zero SNe la near maximum luminosity. The SN2002dc spectrum 
closely resembles that of SN1981B in the galaxy NGC 4536, including the Si II 
absorption at 6150 A, which is the defining feature of SNe la. This is the first 
time this feature has been seen at z 0.5. The spectrum of SN2002dd appears 
intermediate between those of SN1981B and SN1990N in NGC 4639. In par- 
ticular, it shows the deep Ca II trough at 3750 A and numerous other features 
due to Si II and Fe II common to SNe la. 
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Fig. 4. — Observed ACS grism spectrum of the SN2()02dd candidate host 
galaxy (see Fig. 2). The l-cr error region is delineated by the gray area. The 
dark curve is a spectrum of a nearby starburst galaxy from Calzetti et al. (1994) 
that has been redshifted to z=0.95, smoothed to the resolution of the ACS 
grism extraction, and rescaled for comparison purposes. The broad line near 
7200 A in the grism spectrum aligns well with the [O II] A3727 emission in the 
comparison spectrum for this redshift. 
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For SN2002dd, we find good spectral matches with 
SN1995D, SN1990O, SN1981B, SN1972E, SN1990N, and 
SN1994ae, among others. The range of acceptable ages is fairly 
broad, -7 to +9 days past maximum. The best-fit redshift is 
z=0.944± 0.022, with a total range of 0.90-0.97. These re- 
sults rule out the z = 0.79 galaxy l."6 away as a possible host, 
but are consistent with the z = 0.950 redshift of the galaxy 3 "6 
away (our grism spectra for both these galaxies support the 
results from Cohen et al.). We note that a number of spec- 
tral templates, including SN1989B, SN1992a, SN1994m, and 
SN1996x, produced formally acceptable fits, but with zr^O.S. 
In fact, many of the templates have some correlation peak at 
this redshift, but these fits are all of significantly lower quaUty 
than the best fits with z « 0.94. In addition, the photometric data 
for SN2002dd is completely incompatible with zwO.5. While 
the correct redshift is clear in this case, we wish to caution that 
SNe la cross-correlation results in general depend critically on 
the wavelength coverage and the set of spectral templates used. 

Figure 3 shows the calibrated, rest-frame spectra of the su- 
pemovae with a 5-pixel (2-2.5 times the FWHM) extraction 
aperture, and compares them to smoothed spectra of two nearby 
"standard" SNe la at maximum luminosity. Both SN2002dc 
and SN2002dd display many of the standard spectral charac- 
teristics of SNe la, which have been discussed in detail by 
Coil et al. (2000) for SNe la at these redshifts. In particular, 
we note that SN2002dc shows the Si II absorption trough at 
rest A = 6 150 A, which is considered the defining feature of 
SNe la and thus confirms the Type la identification. This is the 
first time the complete feature has been unambiguously seen at 
z « 0.5. It is nearly impossible to observe from the ground at 
these redshifts, being an absorption feature in a faint spectrum 
at A > 9000 A. 

While the SN2002dd spectrum does not extend to the Si II 
A6150 feature, its general characteristics easily rule out all other 
types of supernovae except the rare luminous Type Ic (see Cloc- 
chiatti et al. 2000). However, as discussed by Riess et al. 
(1998), Type Ic supernovae tend to lack the very deep Can 
A3750 absorption with strong recovery to the peak at ~ 4000 A. 
In addition, only SNe la spectra show the notch in this peak due 
to Si II A4130. Clearly, SN2002dd shares these characteristic of 
SNe la. We therefore regard both objects as spectroscopically 
confirmed SNe la having redshifts consistent with the host or 
neighboring galaxies. 

Finally, Figure 4 presents the grism spectrum of the faint 
/ab~25.5 candidate host galaxy identified in Figure 2, and 
compares it to that of the starburst galaxy NGC 3049 from 
Calzetti et al. (1994)'^ The starburst spectrum has been red- 
shifted to z=0.95 and smoothed to the resolution of the grism 
spectrum. The feature at ~7200 A in the candidate host's spec- 
trum matches the [O ii] A3727 fine at this redshift. The contin- 
uum shape is also a reasonable match, although this varies con- 
siderably among starburst galaxies. We suggest that SN2002dd 
lies in the outskirts of this faint galaxy, which in turn is a dwarf 
companion of the larger z = 0.95 galaxy. 

2.3. Photometry 

We measured the fluxes of the two supernovae in each band- 
pass at each epoch using the Vista "psf" routine, which derives 
the scale factors between faint point sources and brighter PSF 
template stars using an implementation of the DoPhot (Shechter 
et al. 1993) "wingy Gaussian" fitting function. Seven different 
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high signal-to-noise (S/N) stars in the field were used as psf 
templates for the F850LP images, and five of these were also 
used as templates for the F775W photometry (the two bright- 
est were just saturated in the first epoch F775W image). The 
aperture flux of the faint source is derived from the aperture 
flux of the high-S/N template times the fitted scale factor. The 
total flux was then derived using the bandpass-dependent aper- 
ture corrections of Sirianni et al. (2003), plus small adjustments 
for the color-dependence of the F850LP aperture correction 
(Gilliland & Riess 2002). We experimented with several differ- 
ent aperture sizes for the PSF stars but always derived consis- 
tent total magnitudes at the 1% level after applying these aper- 
ture corrections. In the end, we chose an aperture of radius 
0''5 as affording the best compromise between systematic and 
random error. Although the F850LP correction becomes very 
sensitive to color at small radii for extremely red objects, the 
adjustment from GilUland & Riess for objects with the colors 
of the SNe and PSF stars is only '--^0.01 mag for this aperture. 
The adopted aperture corrections are then 0.077 mag in F775W 
and 0.11 6 in F850LR 

For SN2002dc, it was necessary to subtract the host galaxy. 
We used the very deep WFPC2/F814W image remapped to the 
scale and orientation of the ACS frames and rescaled the galaxy 
flux to minimize the residuals upon subtraction of the galaxy 
from each of the F775W and F850LP discovery images. The 
flux ratios were determined after smoothing the ACS images 
to the WFPC2 resolution, but the photometric measurements 
were done after subtracting the rescaled galaxy from the origi- 
nal (unsmoothed) ACS images. Figure 1 illustrates our F850LP 
subtraction. Scaling of the flux to the later (shorter) observa- 
tions was done according to the exposure time. We varied the 
derived scale factors by ±0.5% (about twice the uncertainty of 
the flux ratio determinations) and repeated the SN2002dc flux 
measurements. The results differed by less than 0.05%, pri- 
marily because the psf routine always subtracts a locally deter- 
mined sky value. Therefore, the host galaxy subtraction is not 
a significant source of error. 

We empirically determined the magnitude errors and biases 
via Monte Carlo experiments. We constructed a composite PSF 
star for each bandpass and added 50-100 cloned versions of 
this star scaled to the measured magnitudes of each of the two 
supernovae. We then recovered these artificially introduced 
stars, measured their magnitudes with the psf routine described 
above, and compared to the true magnitude. We repeated this 
procedure 10 times at each epoch for both F775W and F850LP. 
The errors determined in this way agreed closely with the errors 
we estimated from first principles based on all sources of noise, 
including the scatter in the results arising from different PSF 
templates. The estimated and measured errors differed on aver- 
age by about 0.01 mag, and never by more than 0.02 mag. We 
found some small biases in the fitted scale factors, ranging be- 
tween zero and 0.06 mag (typically half the measured 1-cr error 
or less), and corrections for these were appUed. We also tested 
our approach against results obtained using DAOPHOT (Stet- 
son 1987), and again found agreement at the 0.01 mag level. 

For die single-epoch NlCMOS/FllOW observation of 
SN2002dd, we took a sUghtly different approach. The small 
field of view encompassed essentially no good PSF comparison 
stars. So, instead, we created ten PSFs for this camera using 
Tiny Tim (Krist & Hook 2001) and added them to each cali- 



brated image after determining the offsets but prior to drizzle 
combination. We then determined the relative scale between 
SN2002dd and these artificial PSFs in the drizzled image using 
botii DAOPHOT and the Vista/DoPhot PSF-fitting approach. 

The ACS instrumental magnitudes were then transformed to 
the AB system using the latest calibrations. This is the same as 
derived from the header photometric keywords for F775W, but 
we anticipate a correction of 0.064 mag to the current (2002 
December) F850LP header calibration (Sirianni et al. 2003). 
The adopted zero points for 1 e"sec"' are then 25.654 (F775W) 
and 24.850 (F850LP)^^. For the hght curve analysis, we require 
Vega magnitudes, and use AB to Vega conversions of -0.401 
(F775W) and -0.569 (F850LP). For the NICMOS/Fl lOW data, 
we calibrated directly to the Vega system using the prelimi- 
nary post-servicing mission zero point for camera 2, which is 
about 55 ±5% more sensitive than for cycle 7 (M.Dickinson 
& M. Rieke, priv. comm.). Because of the still-tentative nature 
of this calibration, we have added a conservative 0.07 mag in 
quadrature to the Fl lOW measurement error to get the final un- 
certainty. 

FinaUy, because the photometric record is most complete 
in F850LP but the first F850LP observation ("epoch 2") oc- 
curred 10 days after the first F775W and G800L observations 
("epoch 1"), we derived a first epoch F850LP magnitude from 
the grism spectrum and F775W magnitude. To do this, we used 
the extracted count spectrum, divided out the G800L response 
function, integrated the result across the filter bandpasses, and 
applied the AB zero points to derive a color directly from the 
grism spectrum, (j-z)gr. Because of aperture and other effects, 
this wiU not exactly correspond to the colors measured from the 
direct images. However, we plotted the epoch 2 F850LP mag- 
nitudes against the epoch 1 F775W-(j-z)gr measurements for 
all isolated, unsaturated point sources having AB mag < 25. 
The result is shown in Figure 5. 

A total of 21 point sources were used to define the Unear re- 
lation in Figure 5, which has a slope of 0.98 and a scatter of 
0.045 mag. The residuals about the relation show no correla- 
tion with (i-z)gr color. Forcing a slope of unity increases the 
scatter by 0.009 mag. Both supernovae lie ~ 0.2 mag above the 
relation, indicating that they each faded by about this amount 
in the 10 day interval. One other relatively faint object (shown 
with an open circle in Figure 5) appears to have faded by a 
similar amount and was excluded from the fit, being a > 4- 
fj outlier, although including it does not significantly change 
the results. We use the fitted relation to "predict" the epoch 1 
F850LP magnitudes of the SNe la, and include a 0.05 mag error 
in quadrature to allow for the scatter about the fit. Our final set 
of photometric measurements transformed to the Vega system 
are presented in Table 1 . 



The 2003 February photometric recalibration of ACS differs from these zero points by -0.013 mag in F775W and -0.006 mag in F850LP. These differences are too 
small to affect any of our conclusions. 
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Fig. 5. — Derivation of the first epoch F850LP SNe la magnitudes. Top: 
Second-epoch F850LP aperture magnitudes are plotted against the difference 
of the first-epoch F775W aperture magnitudes and (i-z)gr color measured by 
integration of the grism spectra. Magnitudes are on the AB system. The soUd 
line was fitted to the 21 soUd points and is used to estimate what the F850LP 
magnitudes of the two supemovae would have been at the first epoch; the off- 
set of the supemovae above the relation indicates that they both faded by about 
0.2 mag in the 10 day interval between epochs. One other object, shown by 
the open circle, appears to have gotten fainter by a similar amount and was not 
used in the fit. Bottom: Magnitude residuals about the fit. The scatter of the 
solid points is 0.045 mag. 

3. LIGHT CURVES AND DISTANCES 

From the multicolor light curves defined by Table 1, we can 
calibrate the SNe la luminosities, and so infer their distances. 
Several different methods exist for doing this, including the 
Amis decline rate parameterization (Phillips 1993), Amis tem- 
plate fitting (Hamuy et al. 1995, 1996), multi-color hght curve 
shape (MLCS) (Riess et al. 1996, 1998), stretch (Perlmutter 
et al. 1997), and best average template (BAT) (Tonry et al. 
2003) methods. Here we employ both the MLCS and BAT 
methods to estimate the distances to SN2002dc and SN2002dd 
and place them on the Hubble diagram. BAT, while less devel- 
oped than MLCS, does not assume a single-parameter family 
of objects, and therefore we can use it to test for consistency 
between two significantly different approaches to measuring 
SNe la distances. 

3.1. K-corrections 

All methods for estimating SNe la distances depend on how 
well the data on the distant program objects and the calibrat- 
ing set of local SNe la can be transformed to a common ref- 
erence frame and photometric system for comparison. This 
can be achieved either by transforming the magnitudes of the 
distant objects to the closest-matching rest-frame bandpass in 



which the calibrating sample have been measured, or equiva- 
lently, through redshifting the photometric record of the cali- 
brators and transforming to the bandpasses in which the dis- 
tant SNe la were measured. In either case, the systematic error 
is minimized if the observed and template bandpasses match 
as closely as possible after shifting to the common reference 
frame. We note that the MLCS and BAT methods used in the 
analysis below take the opposite approaches as far as redshift- 
ing/blueshifting the calibrator/program objects. 

Because the observation and comparison bandpasses will not 
match perfectly even after redshifting to a common frame, it 
is necessary to apply "^-corrections" to convert the observed 
magnitudes to the template bandpasses. The problem of K- 
corrections has been discussed in the context of SNe la by 
Hamuy et al. (1993); Kim et al. (1996); Schmidt et al. (1998); 
and Nugent et al. (2002), among others. We calculate the K- 
corrections in a manner similar to Riess et al. (1998) and Tonry 
et al. (2003), except we use the model spectral energy distribu- 
tions (SEDs) from Nugent et al. (2002). However, we force the 
models to agree with the average SNe la light curve from Jha 
(2002), which is based on extensive empirical multi-color data. 

For each day that the model and template Ught curve data are 
tabulated, we integrate the model SEDs through the U,B,V,R,I 
bandpasses, fit a spline to the differences with respect to the em- 
pirical Ught curve, multiply the result into the SED, and iterate 
to convergence. Once the model SED colors have converged to 
the observed light curve colors for every date, they can be used 
for the /T-correction calculation. In this way, we preserve both 
the detailed spectral features of the models and the general con- 
tinuum shape defined by the empirical light curves. We include 
an allowance for ±0.07 mag error in the /sT-corrections based 
on the scatter derived from individual SNe la near maximum- 
light (B.Schmidt, priv. comm.). This inethod ensures consis- 
tency with previous SNe la distance analyses while allowing us 
to obtain reUable A'-corrections for our rest-frame U -band data. 
We note that if there is significant uncertainty in the redshift, 
this should also be included in the error analysis, but the K- 
correction near maximum is not highly sensitive to the precise 
redshift. Also, in the present case, we beheve that the SNe la 
redshifts are accurate to 0.01 (better for SN2002dc), based on 
a combination of the SNe la and galaxy spectra. Table 1 gives 
the F850LP /T-corrections used in the MLCS analysis. 

3.2. Distances 

We estimate the SNe la distances using both the MLCS 
method described in detail by Riess et al. (1996) and developed 
further by Riess et al. (1998), and the BAT method described 
by Tonry et al. (2003). Both methods use Ubraries of multi- 
band SNe la hght curves. The main difference is that MLCS 
parameterizes the light curve shape as a continuous 1 -parameter 
function of the luminosity excess A, whereas BAT determines 
the SNe la luminosity from a weighted average among the best- 
matching comparison Ught curves and uses the scatter in its er- 
ror estimate. Both methods fit for the time of maximum and 
yield estimates of the extinction based on the color difference 
between model and observations. The distance is then inferred 
from the caUbrated light curve after applying the derived rest- 
frame extinction correction. This should not be confused with 
the extinction correction of Ay = 0.04 mag (Schlegel et al. 1998) 
that we apply first to the observed magnitudes prior to either 
method. 

Table 2 Usts the results of our MLCS and BAT analyses. The 
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columns are: (1) SNela designation; (2) redshift; (3) MLCS 
luminosity parameter A (mag); (4) MLCS estimated distance 
modulus; (5) MLCS fitted V-band extinction (mag); (6) best- 
matching SNe la light curve template from the BAT method; (7) 
BAT estimated distance modulus; (8) BAT fitted V-band extinc- 
tion (mag). The sense of the A parameter from the MLCS fit is 
that SN2002dc is found to be underluminous by 0.1 mag, while 
SN2002dd is overluminous by 0.4 mag (i.e., systematic errors 
of these sizes would be made in assuming that SNe la were per- 
fect standard candles). The light curve fits for SN2002dd were 
done both with and without the F775 data, corresponding to 
the rest-frame U band. The two methods produce very simi- 
lar results in each case, well within the quoted l-a fit errors. 
Although SN2002dc is nearer and brighter, the estimated er- 
rors from both methods are larger for this object because it is 
found to have significantly higher extinction, which adds to the 
distance uncertainty. The extinction is consistent with the pres- 
ence of SN2002dc within the disk of its host galaxy. 




-20 20 40 60 

observed days 



Fig. 6.— Light curves and fitted MLCS model for SN2002dc in the rest- 
frame V and R bands. The model includes: luminosity parameter A, which 
determines the shape of the curves; reddening, which sets the extinction; time 
of maximvmi; and the distance. All points were used in deriving the fit. 
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Fig. 7.— Light curves and fitted MLCS model for SN2002dd. Only the 
F850LP and Fl lOW measurements (rest-frame B and V, respectively) were 
used in deriving this fit. The rest-frame U data (observed F775W) points were 
not used, but are shown in comparison to the U Ught curve of the fitted model. 
The model appears to decline more rapidly than the data in U, which could be 
due to metallicity or other effects that become much more pronounced in the 
ultraviolet, or due to insufficient [/-band template data. 



The light curve fits for SN2002dc and SN2002dd from the 
MLCS method are shown in Figures 6 and 7, respectively. 
Overall, the fits look quite good. For SN2002dc, the observed 
F775W and F850LP observations are transformed to rest-frame 
V and R, which provide the best matches at this redshift. The 
results are basically unchanged if the transformation is instead 
done to B and V, which is an indication of the consistency of 
the A'-corrections. The /?-band fight curve (observed F850LP) 
shows a hint of the late-time ingress of flux at red wavelengths, 
a common characteristic of SNe la. 

For SN2002dd, the F850LP and FllOW bandpasses fi-ans- 
form most readily to rest-frame B and V, which are the most 
commonly used bandpasses for fitting SNe la light curves. The 
empirical library of U -band light curve data is much less com- 
plete. Moreover, the shape of the UV continuum becomes sen- 
sitive to the metalhcity (e.g., Hoflich et al. 1998), which can 
cause systematic errors in the fits. These problems are com- 
pounded by the dust obscuration being greater in U, so the 
peculiar systematic effects of this poorly constrained bandpass 
can unduly drive the fitted extinction values. For these reasons, 
rest-frame UV data have not in the past been used for light curve 
fitting, and our results omitting the F775W data may be more 
reliable; however, the analysis is not overly sensitive to their 
inclusion. Figure 7 shows the MLCS model from the fit to the 
F850LP and FllOW data alone, but compares to the full data 
set. The agreement is reasonably good in U, indicating that 
SN2002dd is not particularly abnormal in the UV, although the 
data appear to decline more slowly than the model, possibly 
due to the effects mentioned above. A better understanding and 
more extensive UV light curve data sets are needed to take fuU 
advantage of optical observations of SNe la at z > 1 . 

3.3. The Hubble Diagram 

Figure 8 places the new SNela distance measurements on 
the differential Hubble diagram of pubhshed SNe la. The dis- 
tances are relative to a fiducial "coasting," or "empty," universe 
(f7n,=OA=0), thus removing from consideration. The pub- 
lished data are binned as in Riess et al. (2001) with the ex- 
treme SN1997ff point represented by its upper limit. The data 
are compared to expectations from various models, including 
Einstein-deSitter (fim = 1), an open universe model (Jim = 0.3, 
f^A = 0), a flat A-dominated model {fl^ = 0.3, f^A = 0.7), and 
some simplistic intergalactic "dust" models (described below). 
We use the MLCS distances for this comparison because the 
published data in the figure are largely from the MLCS method 
or have been homogenized to remove any mean offset from 
MLCS. The data point for SN2002dc would be offset fainter 
by 0.17 mag using flie BAT result, while the SN2002dd point 
would be virtually unchanged. 

SN2002dc, like most individual SNe la at similar redshifts, is 
consistent with any of these models except ilm= 1 without dust. 
SN2002dd at z=0.95 favors the A-dominated or dust models, 
but SN1997ff at z ~ 1 .7 remains the sole indicator of past decel- 
eration. A constrained JIa + = 1 fit to the binned data (open 
symbols) gives Oa = 0.74 ±0.08; adding the two new MLCS 
distances and performing, purely for illustrative purposes, a hy- 
brid binned/unbinned fit results in JIa = 0.73 ± 0.07. If we re- 
place the MLCS distances with our BAT results, the value be- 
comes fiA = 0.74 ± 0.07, or = 0.75 ± 0.07 if the SN2002dd 
distance including the rest J/-band data is used. Thus, our data 
are fuUy consistent with the findings of Riess et al. (1998) and 
Perlmutter et al. (1999) favoring an accelerating universe. 
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Fig. 8. — The differential Hubble diagram of previously published and new SNe la compared to various cosmological models, with and without grey dust obscu- 
ration. SN2002dc and SN2()02dd from this work are shown with solid symbols (MLCS distances and errorbars). The previously published, mainly ground-based 
data have been grouped into six redshift bins as in Riess et al. (2001) and are shown with open circles. The original z<0.1 data come from Hamuy et al. (1996) 
and Riess et al. (1999), while the 0.3 J; z ^ 0.8 data are from Riess et al. (1998) and Perlmutter et al. (1999). The extreme point at z= l .V indicates the upper 
limit for SN1997ff from Riess et al. (2001). The "grey dust" model represented by the long-dashed line has a constant comoving dust density and is similar to 
the one considered by Riess et al., except here we have calculated it for the open r2,n = 0.3, Sl\ =0 universe, rather than for an "empty" one. The "evolving dust" 
models shown by the heavy (Sim = 0.3) and light (fim = 1 0) short-dashed lines have a comoving dust density that builds up with decreasing redshift to compensate 
for dilution from cosmic expansion. The dust models are arbitrarily normalized to go through the data points at z~0.5. See text for details. 



However, Aguirre (1999) has argued that a uniform distri- 
bution of intergalactic dust can explain the faintness of inter- 
mediate redshift SNe la without invoking cosmic acceleration. 
Such dust might be expected to have grey scattering proper- 
ties because only large grains (> 1/im) would likely survive 
the extreme conditions required for ejection into the intergalac- 
tic medium and the high temperatures believed to prevail there. 
While the observed brightness of the z ~ 1.7 SN1997ff argues 
against the simplest, non-evolving, grey dust model and in favor 
of cosmological deceleration beyond z= 1 (Riess et al. 2001), 
Goobar et al. (2002) show that models with an evolving dust 
density can accommodate SN1997ff. 

In Figure 8, the model labeled "grey dust" is very similar to 
the one considered by Riess et al. (2001), with the optical depth 
directly proportional to redshift to approximate a fixed comov- 
ing density of dust, except here we calculate it for the = 0.3 
open model, rather than for = 0. However, the behavior of 
this non-evolving grey dust model is fairly insensitive to fim- 
The normalization of the extinction is adjusted to match the ob- 
served SNe la data near z ~ 0.5. The "evolving dust" models 
shown in the figure use the same cosmology, as well as f2ni= 1, 
but have constant spatial densities of dust, to approximate con- 
tinual injection of dust at a high enough rate to compensate for 
the dilution effect of the expansion; thus the comoving density 
decreases to higher redshift. 

These two types of dust model essentially reproduce the be- 
havior of the Goobar et al. (2002) models A and B, which 
those authors admit are not necessarily plausible in themselves 
but provide "useful limiting cases." The plotted "dust" curves 
can therefore be taken to mark the approximate range of high- 
redshift behavior for decelerating models which invoke grey in- 
tergalactic dust. Of course, perfectly grey dust does not exist in 
reality, and efforts to detect a wavelength dependence of the 
SNe la dimming have not found any (e.g., Riess et al. 2000). 
Moreover, a direct upper limit on the density of intergalactic 



dust derived from the absence of a detectable X-ray scattering 
halo around a single distant quasar is well below that needed to 
explain the dimming of SNe la (Paerels et al. 2002). Therefore, 
cosmic acceleration seems the most natural and compelling ex- 
planation, although only further observations of high-z SNe la 
over a large wavelength range can definitively rule out the chal- 
lenge posed by these dust models. 

3.4. Systematic Effects 

Various systematic effects on high-z SNe la distances such as 
extinction, gravitational lensing, and selection bias have been 
discussed in detail by Riess et al. (1998), Perlmutter et al. 
(1999), and several subsequent works. In our case, the ob- 
jects were not found as part of any well-defined search effort, 
and so a statistical correction for selection bias is not possible. 
However, because they were both found well above the detec- 
tion limit (the statistical error at the discovery epoch being only 
0.01-0.02 mag), this bias would be minimal. We note that pre- 
vious high-z SNe la searches have not used such corrections. 

Two issues are worth noting, however. First, the fitted extinc- 
tion for SN2002dc, although fairly modest in Galactic terms, is 
large compared to most other high-z SNe la previously mea- 
sured. It is higher than for the great majority of objects in both 
the Riess et al. (1998) and Tonry et al. (2003) samples (most 
of which have fitted extinctions consistent with zero), as well 
as most objects from Perlmutter et al. (1999). On the other 
hand, the uncertainty in the extinction is fairly large, ±0.19 
mag from MLCS and ±0.30 mag from BATM. This is the dom- 
inant source of distance uncertainty for SN2002dc (the extinc- 
tion error for SN2002dd is about ±0.16 mag from both meth- 
ods). Still, it may be that seeing effects in ground-based data 
tend to select against SNe la deep inside galaxies that would 
have higher extinction values. Therefore, the quality of highly 
extincted SNe la distances has not been thoroughly tested. Fur- 
ther, there are no external checks of the extinction estimate. 
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which depends critically on the assumption that the redden- 
ing law in the z=0.47 host is the same as that for the interstel- 
lar medium of our Galaxy, but variations are known to occur 
(CardelU et al. 1989; Falco et al. 1999). 

Another effect worth considering is weak gravitational lens- 
ing. Stochastic lensing by large-scale structure will tend to 
demagnify distant SNela because the clustering of galaxies 
causes the typical line of sight to be underdense (e.g., Kan- 
towski et al. 1996; Wambsganss et al. 1997). Both Riess et al. 
(1998) and Perlmutter et al. (1999) concluded that this effect 
was neghgible for SNela in the redshift range 0.5-1.0. How- 
ever, certain SNe la may be significantly magnified if their line 
of sight is near massive foreground objects. In other words, al- 
though the magnification probability function peaks at a value 
less than unity, it is asymmetric with a tail towards higher mag- 
nification. In the case of SN1997ff at z« 1.7, which is pro- 
jected near several bright z = 0.5-0.9 galaxies, the magnifica- 
tion has recently been estimated at -0.34 ± 0.12 mag (Bemtez 
et al. 2002), which would bring it very close to the "empty" 
model in Figure 8. 

We performed the same analysis for the two new SNela as 
done by Bemtez et al. (2002) for SN1997ff. The magnification 
derived for SN2002dc at z « 0.5 is negligible. However, we find 
a magnification of -0.22 ±0.10 mag for SN2002dd, with the 
largest contributions coming from galaxies 381 and 345 in the 
Fernandez-Soto et al. (1999) catalogue (3-486.0 and 3-610.1, 
respectively, from Williams et al. 1996). This would mean that 
SN2002dd is ^ 10% more distant than inferred above. Such a 
high magnification at z« 1 would have a probability of < 0.5%, 
based on the A-CDM galaxy-lensing model of Holz (1998), and 
an even lower probabiUty according to the nearest equivalent 
model from Wambsganss et al. (1997). However, it is only sig- 
nificant at the ~ 2-<T level, and this analysis does not account for 
the demagnification tendency of large-scale structure, which is 
very difficult to estimate for individual objects. Thus, our mag- 
nification estimate is biased high at some level. 

To explore this issue, we repeated the analysis for 1000 ran- 
dom lines of sight in the HDFN and found an average magni- 
fication of -0.05 mag, with 1.8% of them yielding values as 
high as for SN2002dd. If the average is taken as a zero point 
offset, the magnification for SN2002dd becomes -0.17 ±0.10 
mag. Of course, any given estimate may have systematic errors 
from the uncertain mass-to-light ratios of the foreground galax- 
ies and the shapes of their potentials. Magnification corrections 
have not been applied to other SNe la at z ^ 1 , and although we 
also choose not to apply this correction, we note it is probable 
that SN2002dd is magnified and the effect in Figure 8 would be 
to move this data point upward by perhaps ~ 0.2 mag. 

4. DISCUSSION AND SUMMARY 

Our analysis has added two additional high-z points to the 
SNela Hubble diagram, the first that have been found with 
the new ACS instrument on HST. Both objects were spectro- 
scopically confirmed as SNela near maximum-hght by ACS 
grism spectra. SN2002dd at z=0.95 has a higher redshift than 
any other spectroscopically confirmed SNe la with a published 
distance. For instance, SN1997ck at z=0.97 (Garnavich et al. 
1998) had neither spectroscopic confirmation nor useful color 
information; SN1997ff at z«1.7 also lacks spectroscopic con- 



firmation, has an uncertain redshift, and all the photometry was 
transformed to the B-band, so there was no constraint on ex- 
tinction; SN1999fv at z« 1.2 is the most distant SNela with 
pubUshed confirmation (Tonry et al. 1999; Coil et al. 2000), but 
has insufficient light curve data with no color information, and 
thus only a weakly constrained distance (Tonry et al. 2003). 
However, with the already heavy use of ACS by the SNela 
community, we expect that this circumstance will not last for 
very long. 

Taken at face value, our distance for SN2002dd is most con- 
sistent with the A-dominated cosmology favored by the analy- 
ses of the larger samples of ground-based SNe la in the z = 0.3- 
0.9 range by Riess et al. (1998) and Perhnutter et al. (1999). 
Within the context of such flat models dominated by a cosmo- 
logical constant, the transition from acceleration to decelera- 
tion occurs at the "coasting point" redshift Ze = [2f2A/f^m] - 1 
(Turner & Riess 2002). A reUable measurement of this redshift 
would therefore provide a strong constraint on rt\. Although 
the transition occurs at z ~ 0.7 for models with A « 0.7, the red- 
shift of SN2002dd is still too low to yield any clear evidence for 
deceleration; in particular, it remains consistent with the sim- 
plest grey dust model. Moreover, we estimated that SN2002dd 
may be magnified by ^0.2 mag due to the lensing effect of 
foreground galaxies, and this would push it further into agree- 
ment with the range of low-fim dust models. A larger sample of 
similar quality data on SNela reaching to z ^ 1.5 is needed to 
explore fully the systematic effects, definitively test the alterna- 
tives to acceleration, and constrain the expansion history of the 
universe. 

In our view, the most encouraging aspect of the present work 
is the relative ease with which the observations were accom- 
plished. Two supernovae were found serendipitously in a sin- 
gle ACS frame without even the usual search technique of im- 
age subtraction (the more distant one was found because of its 
grism spectrum; later searching by image subtraction turned up 
no additional bright candidates). Both were confirmed from 
their spectra as Type la near maximum-light and were followed 
with a relatively small expense of time, primarily with ACS. 
We have also shown, using an empirical relationship, how ACS 
grism data can be used to derive accurate colors and improve 
light-curve coverage. Moreover, now that NICMOS is again 
functional, it is possible to sample the rest-frame visual light of 
z> 1 SNela, as we have done for SN2002dd. This aU bodes 
well for systematic searches for high-z SNela using ACS; in- 
deed such a program has already shown extraordinary promise 
(Giavalisco et al. 2002; Riess et al. 2002), and we have every 
expectation that ACS will continue to deliver on this promise. 

ACS was developed under NASA contract NAS 5-32864, 
and this research has been supported by NASA grant NAG5- 
7697 and by an equipment grant from Sun Microsystems, Inc. 
The Space Telescope Science Institute is operated by AURA 
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K. Anderson, J. McCann, S. Busching, A. Framarini, S. Bark- 
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ACS project at JHU, and to Mark Dickinson, Marcia Rieke, 
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Table 1 



Photometry of HDFN Type Ia Supernovae 



Epoch 




F775W'' F850LP'' 


FllOW'' 




SN2002dc 


1 


405.4 


22.39 ±0.02 22.13 ±0.06" 




-0.86 


2 


415.2 


22.33 ±0.02 




-0.85 


3 


438.1 


23.55 ±0.03 23.18 ±0.04 




-0.92 


4 


449.4 


24.05 ±0.06 23 .45 ±0.07 




-0.89 


5 


455.2 


23.82 ±0.06 




-0.86 


SN2002dd 


1 


405.4 


23 .64 ±0.03 23. 19 ±0.06" 




-1.30 


2 


415.2 


23.41 ±0.02 




-1.34 


3 


438.1 


24.98 ±0.06 24.32 ±0.08 


23.89±0.10 


-1.45 


4 


449.4 


25.50±0.10 24.84±0.12 




-1.55 


5 


455.2 


24.96 ±0.09 




-1.57 



"JuUan Date-2,452,000 
''Vega magnitude. 

'^A'-correction used in transforming the observed F850LP (z-band) mag- 
nitude to rest-frame R (SN2002dc) or B (SN2002dd). 

"Reconstructed from the F775W magnitude and G800L spectrum at this 
epoch. See text. 



Table 2 



SNeIa Distances 









MLCS 






BAT 




SN 


z 


A 


{m-M)o^ 


Av 


Tempi. 


(ot-AOo" 


Av 


2002dc 


0.475 


-1-0.10 


42.20 ±0.23 


0.18 


1998aq 


42.37 ±0.35 


0.29 


2002dd 


0.950 


-0.40 


44.15 ±0.20 


0.10 


1995al 


44.27 ±0.27 


0.10 


2002dd, no f/?. . 


0.950 


-0.42 


44. 17 ±0.20 


0.11 


1990N 


44. 17 ±0.27 


0.09 



^Absolute distance moduU for a zero point with Hq =65 .0 km s ' Mpc ' . 
''Results omitting the F775W data (rest-frame U band). 



